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We investigate a class of stochastic aggregation processes involving two types of clusters: active 
and passive. The mass distribution is obtained analytically for several aggregation rates. When 
the aggregation rate is constant, we find that the mass distribution of passive clusters decays alge- 
braically. Furthermore, the entire range of acceptable decay exponents is possible. For aggregation 
rates proportional to the cluster masses, we find that gelation is suppressed. In this case, the tail of 
the mass distribution decays exponentially for large masses, and as a power law over an intermediate 
size range. 
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I. INTRODUCTION 

Aggregation, coagulation, or agglomeration are kinetic 
phenomena in which clusters bond irreversibly to form 
clusters of indefinitely growing mass. A detailed model 
of merging of two clusters into a single cluster should in- 
corporate the mass, position, velocity (or diffusion rate), 
and even geometrical characteristics of each cluster, to- 
gether with the precise merging mechanism. Such a de- 
tailed description is well beyond theoretical analysis. The 
most natural approach assumes that the merging of two 
clusters of masses x and y into a cluster of mass x + y 
occurs with a rate kernel K(x,y). The cluster densities 
then evolve according to the Smoluchowski coagulation 
equations jy . This theoretical framework |l| |s| has found 
numerous applications in physical chemistry as well 
as in astronomy |?]] and mathematics ||. 

Characterizing a cluster by a single parameter, its 
mass, may indeed be a drastic oversimplification. For 
instance, two clusters with a similar mass can be sub- 
stantially different. Additionally, the character of a clus- 
ter may be altered every time it undergoes aggregation, 
and for example, it may become less likely to partici- 
pate in further aggregation events. In this paper, we 
investigate situations where active clusters may become 
passive after merging events. A concrete example where 
active and passive clusters coexist are multi-phase coars- 
ening processes in one dimension. Indeed, upon merging, 
domain walls may remain active or become passive de- 
pending on the phase of their neighboring domains [p|-|IT| . 
Another example is polymerizations of linear polymers. 
Here, if the end monomers can be chemically active or in- 
ert, newly produced polymers may become passive. Scal- 
ing properties and numerical studies of several stochastic 
aggregation problems will be presented in jl2| . 

Specifically, we consider a class of models where newly- 
formed clusters remain active with a fixed probability p; 
otherwise, they become passive in the sense that they 
never aggregate again. As the fate of newly-formed clus- 
ters is determined stochastically, we term this process 
stochastic aggregation. Given a fixed probability for ac- 
tive clusters to become passive, all clusters eventually 



become passive, and the system freezes. Our goal is to 
determine the time dependent distributions of active and 
passive clusters, as well as the final mass distribution of 
passive clusters. 

Within the standard framework it has long been 
recognized that the Smoluchowski master equations are 
tractable for three particular kernels: K{x, y) = 1, % + y, 
and xy j3|. Exact results have also been found for lin- 
ear combinations of the above solvable kernels p3| . [l4| . 
and for one nonlinear interpolation between the constant 
and sum kernels [|l5| . As will be shown below, stochastic 
aggregation is solvable for the three classical kernels. 

For instance, in the simplest case of mass independent 
aggregation rates, we find that Pk, the final mass distri- 
bution of passive clusters, is scale free, i.e, Pk ~ k 1 . 
Since 7 = 2/p, the entire range of decay exponents 
2 < 7 < co is realized by tuning < p < 1. Similar 
behavior emerges in a dual fragmentation process where 
newly-formed fragments may become passive Jl6[ . The 
sum and the product kernels behave quite similarly as 
the final mass distribution is suppressed exponentially 
for large masses. Nevertheless, over an intermediate size 
range, an algebraic decay with a fixed decay exponent 
may occur. The time dependent behavior is governed by 
an algebraically growing size scale, and it can be cast as 
either a scaling or a scaling-like form. In general, the 
growth law for this scale is not universal as it depends 
on the parameter p. 

The rest of this paper is organized as follows. In Sec. II, 
we define the model and outline the general framework. 
Exact solutions of the master equations for the cases 
of constant, product, and sum kernels are presented in 
Sees. Ill, IV, and V, respectively. The results are sum- 
marized and discussed in Sec. VI. 



II. THE MODEL 

Consider the initial conditions where all clusters are ac- 
tive and have the same mass, which can be set to unity 
without loss of generality Then, active monomers ag- 
gregate to form dimers, and so on. After an aggregation 
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event, the newborn cluster remains active with proba- 
bility p, or becomes passive with probability q = 1 — p. 
In the latter case, the cluster does not undergo further 
aggregation. 

Let K(i,j) be the rate by which clusters of mass i 
and clusters of mass j aggregate, and let A k (P k ) de- 
note active (passive) clusters of mass k. The stochastic 
aggregation process can be symbolically written as 



pK(i,j) 



-4; 



Ai + Aj r Sii+j, s±i~TS±j ' JTi+j. 
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We use the notations A k (t) and P k (t) to denote the den- 
sity distributions of active and passive clusters at time 
t. The master equations describing the time evolution of 
the system read 



dA k 
dt 



i-\-j—k 

- y 

2 ^ 



■y 

3>1 



K{k,j)A Jy 



(2) 
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The convolution terms are proportional to p and q, re- 
spectively, reflecting the probabilities of remaining active 
and becoming passive after each aggregation event. The 
master equations are subject to monodisperse initial con- 
ditions 



A k (0) = 6 kil , P fe (0) = 0. 



(3) 



One may verify that the overall mass is conserved, i.e, 

£ fc >iM4fe(*) + *M*)] = i- 

We are interested in the temporal behavior of the mass 
distributions of both active and passive clusters, and in 
the final distribution of passive clusters. An important 
feature of the model is that the number densities of ac- 
tive and passive clusters are ultimately related through- 
out the evolution. Indeed, in each aggregation event, the 
total number of active clusters is reduced by 1 with prob- 
ability p, or by 2 with probability q; hence on average, 
it reduces by 1 + q. Simultaneously, each aggregation 
event increases the total number of passive clusters by q 
on average. Thus, the following properly weighted com- 
bination of the number densities A(t) — £fe>i A k (t) an( i 
P(t) = J2k>i Pk{t) is conserved 



qA(t) + (1 + q)P(t) = const. 



(4) 



This second conservation law follows from the rate equa- 
tions as well. However, it holds generally for stochastic 
aggregation processes, even when the master equations 
(||) do not apply. 

Assuming that initially there were no passive clusters, 
P(0) = 0, and setting the initial cluster density to one, 
A(0) = 1, we see that Eq. (||) gives a neat expression for 
the final number density of passive clusters 



P(oo) = 



1 + 9 



(5) 



Additionally, mass conservation implies J2 k kP k (oo) = 1, 
and therefore the average mass of a passive cluster is 
(1 + q)/q. A q^ 1 divergence occurs in the q — > limit, 
corresponding to traditional aggregation. Remarkably 
the final density (or equivalently, the average mass) is 
independent of the details of the model including the ag- 
gregation kernel, the clusters' transport mechanism, or 
the system's dimensionality. 

In what follows, we solve the master equations for 
three different kernels: constant K(i,j) — 2, product 
K(i,j) — ij, and sum K(i,j) = i + j. Although we 
restrict our attention to monodisperse initial conditions, 
our techniques apply to arbitrary initial conditions. Fur- 
thermore, the nature of the solutions extends to compact 
initial distributions. 



III. CONSTANT KERNEL 

The constant kernel is the most widely used one. It 
has been applied to diffusion limited coalescence, for ex- 
ample. In traditional aggregation with K(i,j) = const, 
the typical mass grows linearly with time, and the mass 
distribution is exponential. We shall see below that the 
latter assertion also applies for stochastic aggregation. 

We conveniently set K(i, j) — 2, and the master equa- 
tions read 

^=py AiAj-2AA k , 



i+j=k 



i+j=k 



(6) 



The number densities of active and passive clusters can 
be obtained directly. The reaction proceeds through two 
channels: A + A — > A with rate p and A + A — > P with 
rate q, and therefore, the number densities obey 



dA , . ,n dP „n 



(7) 



Solving these equations subject to the initial conditions 
A(0) = 1 and P(0) = gives 



A(t) 



1 



l + (l + q)f 



P(t) 



qt 



l + (l + g)t 



(8) 



Indeed, the relation (Q) between A and P holds, and the 
final density P(oo) = q/(l+q) is recovered. Additionally, 
the density of active clusters decays algebraically with a 
universal exponent: A ~ 

One can also study the mass density of active clus- 
ters, M(t) — ~}Z k> ikA k (t). This quantity has a nice 
probabilistic interpretation: it equals the survival prob- 
ability of an active monomer, initially present in the 
system. In other words, it is the probability that such 
monomer belongs to an active cluster at time t. From 
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Eq. (jg) we find that the mass density evolves according 
to j- t M = —2qAM. Using A(t) from Eq. (|) and the ini- 
tial condition M(0) = 1, we integrate the rate equation 
to find 



M(t) = [l + (l + g )t]-rfe. 



(9) 



Hence, the decay of the total mass density and the growth 
of the average cluster size exhibit non universal behav- 
ior as they depend on the parameter q: M ~ £- 2 9/(l+9) 
and (k) = M/A ~ t p/{1+q \ respectively. Note also that 
the linear average mass growth, (k) ~ t, is recovered for 
the traditional aggregation process p = 1. Similarly, the 
monomer density A\(t) = [1 + (l + q)t]~ 2 '( 1+q \ obtained 
from jfiAi = — 2AAx, exhibits an algebraic decay with a 
g-dependent exponent. This quantity is the probability 
that a monomer avoids aggregation events up to time t. 

We solve for the mass distribution of active clusters 
first. To this end, we transform the distribution Ak and 
introduce a modified time variable r as follows 



A k = A k exp 



t = p dt\ exp 
Jo 



2 / dt' A{t') 
o 

ti 

2 / dt 2 A(t 2 ) 



(10) 



The modified time variable can be written explicitly, 
t = 1 — [1 + (1 + q)t]~~ . In terms of this time variable 

one has A = (1 - r) 1 ?, and A k = A k (l - t)~ 2 /p. The 
above transformations reduce Eq. (m to 



dA k 
dr 



(11) 



i+j=k 



These equations can be solved by the generating func- 
tions technique. Indeed, A(z, r) = J2 k >i ekz Ak(r) obeys 



8A(z,t) 
dr 



A\z,t). 



(12) 



The solution is given by A(z,t) = A (z)[l — tAo(z)]^ 1 , 
with Aq(z) = A(z,0) — A(z, 0) the initial generating 
function. For the monodisperse initial conditions (^), 



Ao(z) = e z , and we find A k 
mass distribution reads 



Thus, the original 



A fe (r) = (l-r) 



(13) 



i.e., it remains exponential throughout the evolution. Al- 
though we present this solution as a function of r, it is 
an explicit solution as r(t) is known. 

The passive cluster distribution can be obtained by in- 
tegrating the corresponding master equation. Again, it 
is more convenient to use the modified time variable as 



Eq. (g) simplifies to ±P k 



(k - 1)(1 - r) 2 /fT fc - 



Integrating over the modified time gives 



P k (r) = -(k-l) dx(l-x) 2/p x k - 2 . (14) 
P Jo 



The final distribution is obtained by setting r = 1. It 
can be expresses in terms of the Euler Gamma function, 



q r(i + 2/ P )r(fc) 
Pfc(oo) = P r(k + 2/ P ) ' 



(15) 



The large mass tail of the distribution is therefore sup- 
pressed algebraically according to 

P k {oo) ~qp- 1 T(l + 2/p)k- 2 /P. (16) 

Interestingly, given the algebraic decay P k ~ k^ 1 as 
k — > oo, mass conservation restricts the exponent range 
to 7 > 2. In our case 7 = 2/p, and since < p < 1, the 
entire range of acceptable exponents emerges by tuning 
the only control parameter p. This behavior is reminis- 
cent of a related stochastic fragmentation process where 
newly-formed fragments may turn stable [Tq]. 

Despite their different limiting behaviors, both mass 
distributions obey the same scaling form. This reflects 
the nature of the process: the distributions are coupled, 
and the "activity" in the system involves masses of the 
order of the characteristic mass. To obtain the scaling 
forms, we consider the limit k — > 00 and t — > 00, with 
the scaling variable £ = k{l — r) kept finite. In this limit, 
Eqs. (|j"3|) and (O) acquire the following scaling forms 



A k {t) ~ t-—«F(0, Pk(t) ~ t- — G ^) 



with the scaling functions 



(17) 



F(0=exp(-0, G(0=r 2/p r(l + 2/p,0. (18) 

Here, T(a, £) = /"dij;"" 1 e" 1 is the incomplete 
Gamma function. The scaling variable £ = fc/fc* contains 
the characteristic mass k* = (1 - r) _1 = (fc) — tP/O+'J). 
At time t, the mass of passive clusters being produced 
is on the order k* . Masses larger than this scale are ex- 
ponentially rare as G(£) oc F(£) = e~^. On the other 
hand, masses smaller than this scale have already turned 
passive. This is manifested by the small argument diver- 
gence G(£) ~ £~ 2 / p , which leads to the time independent 
algebraic distribution Jl6|). 

In fact, the scaling functions are unique. When the 
initial mass is finite, the scaling exponents and conse- 
quently, the scaling form (]l7]) are dictated by the rate 
equations for the number and mass densities. Substitut- 
ing this scaling form into the rate equations yields 
an integro-differential equation for the scaling function 
dr) F(t])F(£ - 77) + = 0. This equation can be 

solved via the Laplace transform, and as long as the first 
two moments of F are finite, the solution is a simple 
exponential. A similar analysis can also be carried for 
the second scaling function. We conclude that the above 
scaling behavior holds in general for compact initial dis- 
tributions. 
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IV. PRODUCT KERNEL 



t = M' 1 - M p/q . 



(24) 



The aggregation process with a reaction rate propor- 
tional to the mass of both reactants, K(i,j) = ij, has 
been applied to polymerizations of branched polymers 
||, and to random graphs |17|. In this case, the sys- 
tem exhibits a phase transition. Specifically, the mass 
condenses into an infinite cluster. In the following, we 
show that gelation does not occur when there is a finite 
probability for clusters to turn passive. 

In this case, the master equations (U) read 

dAk v sr^ a , , a , , 



dt 

dP k _ q 
dt 2 



i-\-j—k 



J2 ^ AiA ^ 

i-\-j—k 



(19) 



Here as well M[t) = J2k>i^ A k(t) is the mass den- 
sity of active clusters. Let us start with analysis of 
the active cluster distribution. To facilitate the solu- 
tion, we again employ the generating functions technique. 
Specifically, the generating function of the sequence kAk, 
A(z,t) — J2k>i kAk(t)e kz , evolves according to 



dA , . , T ,dA 
= (M - M)—. 



(20) 



It is useful to re-write partial derivatives as Jacobians: 

If = and = W§- Substitutin S these Jaco- 

bians into Eq. (|2(]), and using the identity ^ = gfe^ , 
we arrive at 



dz 

at 



= M - pA. 



(21) 



This equation suggests that one should solve for z(A, t). 
Integrating Eq. ( pl| ) gives a solution up to an arbitrary 

constant F(A): z(A,t) = f*dt'M(t')-ptA+F(A). The 
initial conditions (||) read Ao(z) = e z or alternatively, 
F(A) = In A, thereby leading to 



z(A,t) 



dt' M(i') -ptA + hxA. 



(22) 



Therefore, the problem is now reduced to evaluation 
of the mass density of active clusters M(t). Substituting 
A(z = 0, t) = M (i) into Eq. (||) gives 



dt' M(t') - ptM(t) + lnM(£) = 0. 



(23) 



Differentiating this equation leads to (ptM — l)^|f = 
qAl 2 . Rather than solve for M(t), one can obtain an 
explicit solution for t(M). Indeed, the first order dif- 
ferential equation qM 2 j^j — ptM = —1 can be easily 
solved by a number of techniques, e.g., by variation of 
parameters, to give 



As a check of self-consistency, we confirm that this ex- 
pression agrees with the first order Taylor series for the 
mass M(t) = 1 — at, implied by the master equation 
(|l9|). Equation (|24|) in principle gives the mass den- 
sity M(t). To determine the number density A(t), we 
sum up Eqs. © and find jjjA — --^M 2 . Since we 
have an explicit expression for t(M) rather than M(t), 
we treat A as a function of M as well. The density obeys 

which can be integrated to yield 



dA 



= 1+1 
2 



1 + E M x li 
q 



A(M) = — l + q+ pM x l q 



(25) 



Equations (B4h and (|2^) show that both the number and 
mass densities of active clusters decay similarly in the 
long time limit, 



Ait) 



M(t) ~ t- 1 . 



(26) 



The decay exponent is universal and identical to that of 
the constant kernel model. In that case, however, only 
the number density decayed as t -1 . Eq. ( p6] ) also offers 
insight to the nature of the process as the average mass 
approaches a constant (k) — M/A — > 2/(1 + q). One can 
verify that this quantity approaches unity in the limit 
q — > 1 when the active mass distribution contains only 
monomers. 

To complete the solution, we need to determine A k . 
Note that Eq. (S3) can be re- written as Ae~ ptA — 



exp 



f Q dt'M{t') 

is quite handy here 
the equation ve~ v = 
series v = E„m ^ 
tion yields 



. The Lagrange inversion formula 
18 1 . It states that the solution for 



u can be expressed in terms of the 
-u n . Applying it to the above equa- 



A k (t) 



(kpt) 



fe-i 



k ■ k\ 



■ exp 



-k I dt'M(t') 



(27) 



We again express A k as a function of M rather than 
t. The integral can be read from Eq. (|23|), while the 
time t(M) can be read from Eq. (pMj). Substituting the 
corresponding expressions yields the active cluster mass 
distribution in terms of the variable v = M x / q 



A k (v) = v l 



[kp(l-v)] k 1 kv(1 _ v) 
k-k\ 



(28) 



In the large mass limit, the leading behavior of this mass 
distribution is given by 



A k {u) 



yj2np 2 {l-v) 2 



(29) 



with the decay constant A(i/, q) = qv — [q + ln(l — 
l)] — [ v + m (l — v )\- This constant must be posi- 
tive and this is confirmed by the expansion A(v, q) = 



4 



1 V + Er 



1 (q n + v n ). Therefore, the mass distribu- 



tion decays exponentially with mass. This in particular 
implies that no gelation occurs. However, over an in- 
termediate range k <C A _1 (^, q), the mass distribution 
decays algebraically, Ak ~ £ _1 fc~ 5 / 2 . Furthermore, the 
monomer density equals Ay = M e ~v+P M /q j and in the 
long time limit A\ ~ e~ p t~ x . 

The asymptotic behavior is of special interest. In the 
limit k — > oo and t — > oo with the scaling variable £ = kqv 
kept finite, the mass distribution of active clusters can be 
re-casted into a "scaling-like" form 



(30) 



with A(q) = — ln(l — q). The scaling function is simply 
exponential, F(£) = e~^. The above differs from an ordi- 
nary scaling form in that it is cannot be reduced from two 
variables (fc, t) to one rescaled variable. Furthermore, the 
exponential factor e~ Afc indicates that the mass distribu- 
tion is appreciable over a fixed mass range of the order 
1/A(q). This range is significant only in the (traditional 
aggregation) limit q — > 0. In this case, a universal alge- 
braic decay Ak ~ t^fc -5 / 2 is realized in the intermediate 
mass range k -C q~ 2 . Despite this, it is interesting that 
the distribution is characterized by a scaling-like factor, 
with a scaling variable £ = fc/fc*. Note that the char- 
acteristic mass grows with time, k* ~ q -1 ^ 1 / 9 , although 
the average mass (fc) remains finite. 

Turning to passive clusters, we first note that their 
number density can be obtained using A(M) of Eq. (pj£ 
and the conservation law (El) 



limit. Introducing the variable y = xk, keeping only the 
leading terms, and performing the integration yields 



V 



fcA(g) 



(33) 



with A(g) as in Eq. (p0|). Regardless of g, the tail of the 
distribution decays exponentially. However, in the limit 
q — » 0, a universal algebraic decay Pfc(oo) ~ fc~ 3 / 2 occurs 
in the intermediate regime fc <C q~ 2 . The average mass 
divergence q -1 agrees with Eq. (^). Hence, the behavior 
is fundamentally different than that found for mass in- 
dependent rates as the distribution of masses of passive 
clusters is characterized by a finite scale. 

The passive mass distribution d32| ) can be re-cast into 
a "scaling-like" form. Using the same scaling variable 
£ = kqv, one finds 



P k (t) ~ P fc (oo) 



r(«) ' 



(34) 



This form involves the incomplete Gamma function, as 
was the case for the function G(£) in the case of constant 
kernel (see Eq. (jig)). Despite their very different nature, 
the mass distributions of passive clusters in both constant 
and product kernel models have a surprising similarity. 
Indeed, the former distribution (14) can be written in the 
form (pj) with q replaced by I + 2 /p. 



V. SUM KERNEL 



P{M) 



q qM 
l + q ~ ~ 



P 



-M 1 ' q 



(31) 



Alternatively, this density can be obtained by integrating 
j- t P = \qM 2 . While the final number density (H) is uni- 
versal, the mass distribution of passive clusters is model 
dependent. In the previous case of the constant kernel, 
we evaluated the convolution term in the master equa- 
tion explicitly. In the present case, it is more convenient 
to circumvent this sum by noting that it appears in both 
evolution equations. The evolution equation for Pk(t) 
can therefore be simplified to j- t P k {t) = ^[-^A k + kA k M] 
for k > 2. Again, we present the mass distribution in 
terms of v — M x l q 



P 



A k (u)+ / dx 



q + px [kp(l — x)] 



fe-1 



-kp(l—x) 



(32) 



The final distribution P k (oo) can be written explicitly 

using f*dxx l - 1 (l-x) k - 1 a l3x = B(l,k) iF x (/,/ + k, 
with B(l,k) = T(l)r(k)/T(l + k) the Euler Beta func- 
tion and ii 7 i(a,6, x) the confluent hypergeometric func- 
tion jl9| . Additionally, it is possible to evaluate the large 
mass asymptotic behavior as the above integral is dom- 
inated by a region of the order 0(k^ 1 ) near the upper 



The aggregation process with a linear reaction rate, 
K(i,j) = i + j, has been studied because of its math- 
ematical simplicity as well as its relevance to polymer 
formation |2Cfl and to shear-flow coagulation Bj. The 
sum kernel solution is described concisely since it involves 
techniques similar to those used above. The mass distri- 
butions evolve according to 



d -^ = p - k y AiAj 

dt 2 ^ 3 



A k (kA + M), 



i+j=k 



dPk 
dt 



— - k A t Aj 

i+j—k 



(35) 



We introduce the modified densities A k and the modified 
time variable r as follows 



A k = A k exp 



k / dt'A(t')+ / dt' M(t') 



P 



dti exp 



dt 2 M(t 2 ) 



These transformations reduce Eq. 
dA k 



']) into 



dr 



k AiAj 



(36) 



(37) 



i+j=k 
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kz 



The generating function A(z,t) = ^fe>i -4fc( r ) e 
evolves according to §^ = A ^ , which in turn can be re- 
written as |p = —A Thus, for the monodisperse initial 



conditions, we get 



In A - tA. 



(38) 



Again, the crucial part of the solution is the evaluation 
of the number and mass densities A(t) and M(t). These 
are related via the rate equation j^A = —(1 + q)AM, as 
follows from Eqs. (||^). Integrating over time we obtain 



A 1 +i — exp 



dt'M(t') 



(39) 



Inserting this equality into the definition of the gener- 
ating function and using Eq. (|36|), yields A = A~ 
for z = — J Q dt'A(t'). Substituting this equality into 
Eq. (p8|) gives 



lnv4 1 +« - tA^+i 



dt' A(t'). 



Since r = pJ Q dt'A~(t') can be expressed through 
A(t), we conclude that Eq. ( f40| ) is a closed equation for 
A(t). It appears impossible, however, to find an explicit 
solution for A(t) or for A(t). Therefore, we proceed to 
determine t{A). Differentiating Eq. (40) with respect to 
A yields 



1 



1 



,i , - 1 \ q , a dr 
A -tA t+5 I + - A~ —— = 0. 

/ p dA 



(41) 



In deriving Eq. (|41j), we have applied the chain rule 

S-Ilfi aild ^ P J*dt>A^(t>) to evaluate the 
derivative of the right-hand side of Eq. (|0|) . 

We now multiply Eq. (|d]) by A~ ~ and then integrate 
to find an explicit expression for t(A) 



p A~ 



A~ 



(42) 



One can verify that this expression is consistent with the 
expected short time behaviors A = 1— (1 + q)t and r = pt 
implied by Eqs. (35) and (|36|). In the long time limit, the 
above expression yields the leading asymptotic behavior 
of the density and the mass 



A(t) 



1 + 3 



■ r 



M(t) 



1 



1 



r 



(43) 



Hence, these densities decay with the same universal ex- 
ponent as in the product kernel model. Again, the av- 
erage active mass remains finite (k) = M/A — > 1/q for 
t — > oo. One slight difference with the product kernel is 
that the average mass diverges in the limit q — > 0. 

To determine Ah, we again employ the Lagrange inver- 



sion formula. From the exponentiated form of Eq. (38), 
Ae~ rA = e z , we immediately find Ak — (kr) k ~ 1 /k , or 



A k (t) = (krf- 1 exp{- p Q dt'[kA{t') + M{t'))}/k\. Using 
Eqs. ( |40| ) and (|42]), the mass distribution of active clus- 
ters is obtained in terms of the variable v — A 1 ^ 



ifc-1 



kl 



(44) 



Although the variable v is quite different, the similar- 
ity with the corresponding product kernel expression 
( |28| ) is striking. The monomer density in this case is 
A x = Aexp[-p(l - A 1 ^ 1 ^)}. Therefore, it follows the 
same t~ l asymptotic decay as in the product kernel case. 

The appropriate large mass and large time limit is 
again taken in such a way, that the scaling variable 
£ = kqv is kept finite. The corresponding characteristic 
size grows according to k* ~ t 1 ^ 1+q \ The mass distri- 
bution of active clusters can again be re-casted into the 
following scaling-like form 



A k {t) 



(45) 



(40) with A(g) = — q — ln(l — q), and with the scaling func- 



tion F(£) 



The numerous similarities between 



corresponding expressions in the product and sum ker- 
nel solutions reflect the similar underlying mathematical 
structure of both models. For example, the distributions 
generally decay exponentially. Furthermore, in the limit 
q — > 0, there is an intermediate regime k <C q 2 , where 
an algebraic decay Ak(t) ~ fc -3 / 2 occurs. 

We analyzed the mass distribution of passive clusters 
along the same lines as for the product kernel model. For 
example, the final distribution is given by 



Pfc(oo) 



P Jo 



dxx q [(kq-\- 1) + kpx] 



[**>(! -a )] 
fc! 



fc-i 



-kp(l — x) 



(46) 



which may be rewritten explicitly through the beta func- 
tion and the confluent hypergeometric function. Inter- 
estingly, the leading large mass behavior of the final dis- 
tribution is identical to the one obtained for the product 
kernel case (|33|). Additionally, the time dependent be- 
havior can be rewritten using the form ( |34| ) with the ar- 
gument of the Gamma function now being 1 + q (instead 
ofq). 



VI. DISCUSSION 

We have analyzed a class of stochastic aggregation pro- 
cesses. We have solved for the three classical reaction 
kernels: 1, i +j, and ij. Our methods are also applicable 
to the linear combination K(i,j) = a + Hi + j) + cij, 
with a, 6, and c non-negative constants [|l3| , p~4| , and to 
the reaction kernel K(i, j) = 2 — r l — r J , where < r < 1 

We have observed that stochastic aggregation exhibits 
a variety of universal features. For instance, the final 



G 



number density of passive clusters is identical in all three 
models. Also, the number density of active clusters de- 
cays as t _1 . The former feature has been explained by a 
simple argument; remarkably it remains valid even be- 
yond the scope of the rate equations. It would be inter- 
esting to understand the latter feature within a coher- 
ent model-independent framework. For the constant and 
sum kernels, we have seen that both the number den- 
sity, the monomer density and the mass density decay 
according to t" 1 . Furthermore, the results suggest that 
the following scaling-likc form 



A k (t) -put^Fit), 



(47) 



applies when the kernel grows indefinitely with the mass. 
For the solvable sum and product kernels, pk decays ex- 
ponentially in the large mass limit. It will be interest- 
ing to study whether this exponential behavior applies 
for more general kernels, e.g., for homogeneous kernels, 
K(ai, aj) = a x K(i, j), with a positive homogeneity index 
A > 0. 

The simplest constant kernel model is special, and it 
exhibits the most interesting behavior. We have shown 
that in this case, the final mass distribution of passive 
clusters is algebraic. The decay exponent is governed by 
the probability p, and the entire range of decay expo- 
nents, consistent with mass conservation is possible. The 
time dependent behavior follows ordinary scaling, and 
the corresponding mass scale grows algebraically with a 
nonuniversal (p-dependent) exponent. 

Very similar behaviors were found for the product and 
sum kernels. In general, the mass distributions decay 
exponentially, and the average mass approaches a finite 
value. In the limit q — > 0, corresponding to the origi- 
nal aggregation process, an algebraic mass decay occurs 
in the intermediate mass range fc <C q~ 2 . The time de- 
pendent behavior does not follow ordinary scaling, as the 
distributions are dominated by time independent expo- 
nential factors. Nevertheless, a corrective factor which 
is governed by an algebraically growing scale was found 
to be relevant. Remarkably, although the models exhibit 
different temporal behaviors, they possess the same large 
mass behavior of the final mass distribution. 

The above stochastic aggregation provides an inter- 
polation between aggregation (p = 1) and annihilation 
(p = 0) processes. Therefore, it may be used as a tool 
for studying either problems. For example, the Smolu- 
chowski aggregation equations may be applied to low di- 



mensional processes using effective reaction rates. Hence, 
the above solution of the master equations may be rele- 
vant to reaction-diffusion problems. 
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